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Abstract 

We consider multi-particle contributions to nucleon two- and three-point functions from the 
perspective of chiral dynamics. Lattice nucleon interpolating operators, which have definite chiral 
transformation properties, can be mapped into chiral perturbation theory. Using the most common 
of such operators, we determine pion-nucleon and pion-delta couplings to nucleon two- and three- 
point correlation functions at leading order in the low-energy expansion. The couplings of pions to 
nucleons and deltas in two-point functions are consistent with simple phase-space considerations, 
in accordance with the Lehmann spectral representation. An argument based on available phase 
space on a torus is utilized to derive the scaling of multiple-pion couplings. While multi-pion 
states are indeed suppressed, this suppression scales differently with particle number compared to 
that in infinite volume. For nucleon three-point correlation functions, we investigate the axial- 
vector current at vanishing momentum transfer. The effect of pion-nucleon and pion-delta states 
on the extraction of the nucleon axial charge is assessed. We show that couplings to finite volume 
multi-particle states could potentially lead to overestimation of the axial charge. Hence pion- 
nucleon excited states cannot explain the trend seen in lattice QCD calculations of the nucleon 
axial charge. 
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I. INTRODUCTION 


Dramatic progress continues to be made in strong interaction physics by solving QCD 
numerically utilizing Euclidean space-time lattices [l|. Successful determination of the spec¬ 
trum of low-lying hadrons from lattice QCD has been a major triumph [^, and many recent 
calculations are performed at or near the physical pion mass. Small corrections due to strong 
and electromagnetic isospin breaking are becoming relevant for precision determination of 
some quantities, and the recent ab initio determination of the proton-neutron mass split¬ 
ting ( 3 I, for example, shows that lattice QCD techniques are maturing to the point required 
for precision low-energy QCD. Further benchmarks are needed for the nucleon, however, 
before one can reliably use the nucleon as a QCD laboratory for new physics. In this re¬ 
spect, determination of nucleon structure from the lattice remains a prime goal. For a recent 
overview of baryon structure from lattice QCD, see [^. 

The axial charge of the nucleon is one such benchmark quantity, and one for which lattice 
methods have expended considerable effort to determine. This isovector quantity does not 
suffer the complication of quark-disconnected contributions, and so arguably represents the 
simplest aspect of nucleon structure to compute. For pion masses at the physical point, 
however, nucleon correlation functions suffer from a well-known signal-to-noise problem j^, 
which is only exacerbated in the calculation of nucleon three-point functions. For many 
years, lattice computations of the axial charge have been subject to underestimation, with 
a variety of systematic sources investigated for this effect, such as: hnite volume j6l-ll0|. 

1^ 1 and even thermal effects [l^. Not all collaborations. 


excited-state contamination 11 


however, have obtained consistent results, for example, hnite volume effects could not be 
substantiated 151, excited-state effects were not seen in the studies [l^, 1^, and thermal 
effects could not be reproduced [T^. The current status of such lattice QCD computations 


is thoroughly reviewed in 

We investigate the contamination of nucleon two- and three-point functions due to excited 
states. Pion-nucleon and pion-delta excited states are considered within the framework of 


chiral perturbation theory. Some time ago [2^, we determined chiral corrections to the time 


dependence of nucleon correlation functions in inhnite volume by treating the heavy-nucleon 
operator of chiral perturbation theory as a model for a well-optimized lattice interpolating 
held. In fact, the correspondence can be made precise provided the smearing radius of the 
lattice operator is smaller than the pion Compton wavelength, see 12ll. |22| for discussion of 
this and related points. This separation of scales was utilized in |22| to determine model 
independent pion-nucleon contributions to nucleon two-point functions. Chiral perturbation 
theory has also been used to obtain the multi-pion excited-state contributions to the axial- 
axial and vector-vector current correlation functions 23|. In this work, we determine the 


contamination in the nucleon two-point function, and three-point function of the axial-vector 
current by treating chiral dynamics in hnite volume. At leading-order in the low-energy 
expansion, the coupling of pion-nucleon and pion-delta states to these correlation functions 
is predicted on account of chiral symmetry. Such couplings to the two-point function are 
found to be smallj^ and those for the axial current three-point function are shown to favor 
overestimation of the nucleon axial charge. As a result, we conclude that pion-nucleon 
contamination cannot be the source of underestimation of the axial charge. 


This conclusion is also drawn in 2^ . In that work, pion-nucleon contributions to nucleon two-point func¬ 
tions are determined using relativistic baryon chiral perturbation theory. As we utilize a non-relativistic 
expansion throughout, vriTr/MN 1, our results can only be compared in the non-relativistic limit, and 
in this limit all expressions indeed agree. 
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The organization of onr presentation is as follows. First in Sec. m we detail the compnta- 
tion of the nncleon two-point fnnction in chiral pertnrbation theory. One must understand 
how to map lattice QCD interpolating operators for the nucleon into the effective theory, 
and this is done to leading order for the standard (|,0) © (0, |) nucleon operator using 
heavy-nucleon chiral perturbation theory, Sec. Ill A1 Results for pion-nucleon and pion-delta 
spectral weights are obtained on a torus in Sec. IIIBl From the effective mass of the two- 
point function calculated using chiral perturbation theory, we observe that the dominant 
pion-nucleon contamination is quite small, leading to at most a few percent overestimation 
of the nucleon mass. An analysis of the hnite-volume couplings to multi-pion states is car¬ 
ried out using simple A^-body phase-space considerations in Sec. Ill Cl On a torus, we hnd 
the vanishing of spectral weights at threshold, however, not as sharply as in inhnite volume. 
There is one exception, the case of a two-body pion-nucleon state in an s-wave, for which the 
corresponding spectral weight does not vanish at threshold. Away from threshold, one can 
still argue for suppression of all multi-pion nucleon states provided S> 1. Next we in¬ 
vestigate the nucleon three-point function of the axial-vector current in Sec. IIIII The general 
form of the axial-current matrix element is written out in Sec. IIII Al and expanded treating 
the excited-state couplings as small, which is the case in chiral perturbation theory. That 
calculation is carried out in Sec. IIII Bl where contributions from all possible intermediate 
states arising at one-loop order are enumerated. Chiral contamination in the axial-vector 
three-point function is shown in Sec. IIII Cl and drives the correlation function upward, thus 
potentially leading one to overestimate the nucleon axial charge, qa- A summary of findings 
in Sec. nYi concludes our work. 


II. NUCLEON OPERATOR AND TWO-POINT FUNCTION 

We begin by obtaining the spectral representation of the nucleon two-point correlation 
function on a torus using chiral dynamics. To accomplish this, we map the standard lattice 
interpolating field into the corresponding nucleon operator of heavy-baryon chiral perturba¬ 
tion theory. Couplings of pion-nucleon and pion-delta states on a torus are accounted for 
at leading order in the low-energy expansion. We argue for the suppression of multiple-pion 
states using hnite volume phase-space considerations. 


A. Nucleon Operator 

To study the nucleon in lattice QCD, one uses a quark-level interpolating field with 
nucleon quantum numbers, which we label On{x,t). To determine the nucleon mass, for 
example, one computes the two-point correlation function 

G(r) = 5^(0|0„(£,T)0],(0,0)|0). (1) 

X 

Here, for simplicity, we have chosen the nucleon source location at the origin, and projected 
the correlator onto vanishing spatial momentum by summing over all lattice sites, x. In the 
limit of long Euclidean time separation, r, which we have assumed is positive, the two-point 
correlation function has the expected behavior 

G{t) = \Zo\^e-^^^ + ■ ■ ■ , (2) 
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where Mat is the nucleon mass, and the ■ • • represents exponentially suppressed terms. The 
coefficient of the leading exponential, is physically the probability that the interpo¬ 

lating operator Ojv creates the nucleon amidst the tower of excited states, which all possess 
nucleon quantum numbers. In this work, we address contributions from excited states with 
pion-nucleon and pion-delta content. Although we have written the operator 0]y(x, r) with 
a notation that suggests locality, the nucleon interpolating operator need not be local. Gen¬ 
erally one employs some form of gauge-invariant smearing of local operators in order to 
maximize overlap with the ground-state nucleon. As explicated in |22|, the coupling of 
pion-nucleon states to the smeared operator On is hxed by chiral dynamics provided the 
smearing radius remains much smaller than the pion Compton wavelength. 

Written in terms of the quark isodoublet held, q = 
interpolating operator has the formll 

ONi = Qi {q^Cr'^lbq) , (3) 


the standard local nucleon 


where r“ are the isospin matrices. Under a chiral transformation of the form, (L, R) G 
SU{2)l ® SU{2)ji, the quark doublet transforms as the direct sum of two irreducible rep¬ 
resentations, namely (|,0) © (0,|). Given the structure of On in Eq. ([3]), the nucleon 
interpolating operator has precisely the same chiral transformation as the quark isodoublet 
held. Our goal is to hnd an ehective held theory operator for the nucleon that shares this 
transformation. In the context of the ehective theory, such an operator will be local, and 
we necessarily assume the smearing radius for the lattice interpolator is much smaller than 
1/m^. 

In chiral perturbation theory, pion helds are the Goldstone modes emerging from sponta¬ 
neous chiral symmetry breaking, SU{2)l © SU{2)ji —)■ SU{2)y. These modes are embedded 
in the coset held S = exp(2i0 • r//), which transforms as S —)■ under chiral trans¬ 

formations. In our conventions, / = 130 MeV, and the pion helds appear explicitly in the 
matrix </> = 0 • r as 


1 .rrO 


_ { 

TT 


TT 




( 4 ) 


The universal low-energy dynamics of pions is obtained by constructing the most general 
chirally invariant Lagrangian density. This non-renormalizable theory is made tractable 
by utilizing a power counting that treats pion momenta k as small compared to the chiral 
symmetry breaking scale, A^ = 2\/2Trf. Explicit chiral symmetry breaking introduced by 
the quark mass, rrig, can also be included in this power counting. Expanding the coset held 
according to this power counting, S = 1 + 2i(j)/f + • • •, we can see in the ehective theory 
that pions are GauBian huctuations about the chirally asymmetric vacuum. 

The nucleon can be included in chiral perturbation theory as an external source of isospin 
which remains massive in the chiral limit. To maintain a low-energy expansion, one must take 
care of the fact that nucleon momenta p can be on the order of the chiral symmetry brea king 
scale, because Mn ~ A^, where Mn is the nucleon mass. The solution to this problem [26, 27 
is to obtain a non-relativistic formulation by decomposing nucleon momentum modes as 


= Mnv^ + 


( 5 ) 


^ The chiral structure of nucleon interpolating fields has been fully classified in 2^ . This classification has 
been used in the framework of covariant baryon chiral perturbation theory to compute chiral corrections 
to moments of the nucleon distribution amplitude 
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and integrating ont anti-nncleon components of the spinor, which are ~ 2Mi^ away in energy 
from the nncleon components. The resnlting heavy nncleon held, has residnal momentnm 
k, which satishes the power-connting reqnirement k -C A^. 

The form of interactions of nncleons with pions, for example, can be deduced from the 
fact that the nucleon transforms as a doublet under the unbroken SU{2)y symmetry. The 


chiral transformation properties of the nucleon held in the chiral limit are unknown, see [28 


for a nice discussion and conjecture. The inability to resolve soft pions and nucleons in the 
chiral limit, however, gives one freedom to dehne inhnitely many operators which share the 
same SU{2)y transformation, see 29| for a pedagogic overview. The conventional choice is 
to utilize the held 


e = (6) 

which has the chiral transformation ^ , where the unitary matrix t/ is a 

function of L, R, and ^(x). The advantage of this choice is that a nucleon operator living 
in any chiral multiplet can be dressed with pions using the held ^ so that it transforms as 
N.U ^ UNy under SU(2)SU(2)r. In this way, one can construct pion-nucleon interactions 
consistent with the pattern of explicit and spontaneous chiral symmetry breaking without 
knowledge of which chiral multiplet the nucleon belongs. 

In the case of lattice QCD, nucleon interpolating operators, such as that in Eq. ([3]), belong 
to certain chiral multiplets. Given the ehective held theory operator Ny, we must accordingly 
dress it with pions, so that it transforms as (i, 0)©(0, i) under chiral transformations. That 
is we seek an operator N, which can be written in the form Nr + Nr, where Nr —)■ RNr 
and Nr —)■ LNr under chiral transformations. In the low-energy expansion, we have 


Nr = ^^Ny + • • ■ 

NR = ^^Ny + ---, (7) 


where the • • • consist of operators having higher mass dimension. These operators appear 
with unhxed coefficients in the perspective of the ehective theory, and give rise to ehects 
beyond the order to which we are working. The normalization factors of | account for the 
fact that the heavy nucleon held satishes Vr^r'^n = ^Ny + 0{1 /Mn), where Tat is the 
fully covariant nucleon spinor, and the chiral projectors are dehned to be Vr^r = |(1 Tys)!! 
Combining these right- and left-handed doublets into an even parity nucleon operator, we 
see 

Ojr-^N = Zo {Nr + Nr) = ^Zo{^^ + ONy + ---. ( 8 ) 

The mapping of On into N, of course, involves an operator-dependent coefficient Zo that is 
undetermined in the ehective theory, and depends, for example, on the quark-level smearing. 
Our results for pion-nucleon and pion-delta couplings will always be ratios to the nucleon 
coupling, for which \Zci\‘^ cancels out. That this is true can already be seen from expanding 
the nucleon operator N in powers of the pion held 

N = Zo(i--N)n„ + ---. (9) 


^ In terms of the fully covariant spinor 'Pr, we have the right- and left-handed fields Nr = 
and Nr = These lead to a covariant nucleon interpolating operator in the effective field theory 

Nr + Nr = i (^'1' -I- AT -I- 5 which is that found in 221. The 75 -term leads to 0(1/Mr) 

effects in the non-relativistic limit, and these have been dropped in our work. 
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FIG. 1. Additional diagrams contributing to the nucleon two-point function beyond the sunset 
diagrams computed in 20|]. The squares denote the nucleon source and sink locations. The dashed 
lines represent pions, while the solid lines represent nucleons. 


While the overall coefficient of the operator is not hxed, the relative conpling to two pions 
is hxed on acconnt of parity and chiral symmetry. 


B. Nucleon Two-Point Function 


For the nucleon, we have the following Lehmann representation for the inhnite volume 
two-point function 

G{t) = \Zo\^e-^^^ 

with r > 0 assumed for the sink location. This is the large mass limit of the ordinary 
Lehmann representation, with the exception that we work in Euclidean space. The full spec¬ 
tral function, p{E) = S{E) +p{E), includes the nucleon pole contribution to the two-point 
function. The residual spectral function, 'p{E), hence encodes the excited-state contamina¬ 
tion. The long-time limit produces iust the ground-state nucleon, with corrections being 
exponentially suppressed by the gapO 

Using one-loop heavy nucleon chiral perturbation theory, the spectral function for the 
nucleon operator Ny has been obtained previously [1^ by evaluating the sunset diagrams 
shown in Fig. 1 of that work. This result can be written in terms of separate pion-nucleon 
and pion-delta contributions to the residual spectral function 'p{E), namely 


1 + 


dE p{E)e 


-Et 


( 10 ) 


G(t) = \z. 


2 —Mnt 


o e 


1 + 


dEp^N{E)e 


-Et 


+ 


dEp^^{E)e 


-Et 


IrriTy+A 


( 11 ) 


where the pion-nucleon and pion-delta spectral weights are given by 

PnN{ ) ^2 

’ 3(47r/)2 E2 


( 12 ) 


respectively. The characteristic cusp at threshold is universal for a two-body state in a 
p-wave. The gap is m.^ for the pion-nucleon continuum, and -|- A for the pion-delta 
continuum, with A = Ma — Mat as the nucleon-delta mass splitting. 

To derive the spectral representation for the (|, 0) © (0,1) nucleon operator N in Eq. (|8]), 
we must include the contribution from pion tadpole diagrams necessitated by the low-energy 
expansion of the effective field theory operator in Eq. (|9]). Evaluating the diagrams shown 


^ An important distinction is that in Minkowski space, excited-state contributions are only power-law 
suppressed. The leading corrections arise from the accumulation of continuum states near threshold, 
where the spectral function has a cusp. 
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in Fig. [T]in position space, however, we find there is no contribution to the spectral weight 
because the tadpole topology produces only a renormalization of the nucleon wave-function. 
When properly renormalized, this contribution will lead to a correction to the multi-particle 
couplings beyond the order we are working. That the tadpole topology does not contribute 
to the spectral weight can additionally be seen in momentum space. Indeed in Minkowski 
space, there is no imaginary part generated, and hence no multi-particle cut. 

On a periodic spatial torus, the available multi-particle states are discrete. Consequently 
the expected form of the finite volume Lehmann representation is 


G(r) 


\Zo?e-^^" 




\2 —EnT 


(13) 


where u is a collective label for the higher-lying states, and the energies, En are measured 
relative to the nucleon mass. To determine the weights we must revisit the computa¬ 

tion on a torus, which we take to have length L in each of the three spatial directions. We 
keep the length of the time direction infinite throughout. 

For the pion-nucleon contribution, we can return to an intermediate step of the compu¬ 
tation by noting that the energy, is given by E'^ = k'^ + rri^. Hence, we can revert the 
energy integration to the form of a three-momentum integral 



and hence have the replacement 



L3 




(15) 


where the sum is over periodic momentum modes, k = ^n, indexed by the mode number n. 
The energy factor in the denominator is obviously a manifestation of the Lorentz invariant 
measure. Because of the factor \k\ appearing in the denominator, the two-body phase space 
near threshold will thus be removed from the spectral functions above when converting from 
a continuum to a discrete set of states. We will investigate the case of an iV-body final 
state in finite volume below. For the pion-delta contributions, the energy is given in terms 
of the pion three-momentum by the relation {E — A)"^ = k"^ + m^. These facts can be used 
in conjunction with Eq. flT2l) to determine the finite volume couplings to pion-nucleon, and 
pion-delta states. Using the mode numbers as labels, we hav^ 



2 


^9l P 

[^2 _|_ | 722 j 3/2 

24v _ 

\Jk'^ + ml + A 




k2 _|_ ^2 


(16) 


® Our result for which is non-relativistic, agrees with the first term in the expansion of the 

relativistic result derived in 


22 |. 
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FIG. 2. Effective mass of the nucleon two-point function computed with chiral perturbation 
theory, Eq. ()17p . Pion-nucleon excited-state contamination is shown at the few-percent level for 
typical values of lattice parameters. Eor times such that r > 5, the correlator is saturated by 
pion-nucleon states with momentum mode numbers satisfying < 4. 


As we anticipated, these spectral weights vanish with nnlike their continnnm connterparts, 
which vanish as \k\^. Moreover the dependence is now analytic near threshold. 

We can now write the resnlting nncleon two-point fnnction in hnite volnme as 


G(r) 


Zq 


Sg-Mjvr 



(17) 


In the above expression, Efi is the pion energy. As a resnlt, the energy appearing in the 
second exponential, Efi + Mjsi + A = Ea + Ma, is the total energy of the pion-delta inter¬ 
mediate state in the non-relativistic limit. Dne to the additional exponential snppression, 
6 “^’’, contribntions from delta intermediate states are qnite small. The pion-delta spectral 
weights, however, are needed below to calculate the three-point function, for which delta 
contributions are non-negligible. To assess the size of pion-nucleon contamination in the 
two-point function, we form the effective mass of the chiral perturbation theory correlator 

Meff(r) = - log (18) 

G(r) 

where a is the lattice spacing. The effective mass is shown in Fig. |2l where we have employed 
the lattice parameters: a = 0.1 fm, and L = 48 a. We have chosen to plot results at the 
physical value of the pion mass, using the couplings qa = 1-25 and s^atv = 1-5, along with the 
nucleon-delta mass splitting A = 0.29 GeV. The pion-delta contributions, however, cannot 
be discerned in the plot. Pion-nucleon contamination is seen to be at the level of a few 
percent. The smallness of such contributions owes to three facts. 

1. Pion-nucleon interactions are relatively weak, which is exhibited by the factor (/L)“^ 
in the spectral weights. 

2. Lowering the pion mass at hxed lattice size reduces the energy gap, however, it also 
brings one closer to threshold, for which the spectral weights vanish. 

3. Raising the lattice volume at hxed pion mass also reduces the energy gap, but similarly 
brings one closer to threshold. 
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Despite the fact that spectral weights vanish less rapidly near threshold compared to the 
case of a continnnm of states, the snppression at threshold in hnite volnme is enongh to 
ensure pion-nucleon contamination is likely under control for compact lattice interpolating 
operators. 


C. Multi-Particle Contributions 


For a continuum of states, one can derive the behavior of the spectral function near an 


A^-particle threshold [3^. Let the energy of the A^-particle state be F’tk = where 

the Mi are masses of the individual particles. The spectral function then has the behavioiEI 


Pn{E) oc \/E — E' 


-3N-5 


Th 


(19) 


for E > F^Th- From our discussion above, we know this characteristic factor arises from the 
integrals over momentum space, and so must be reconsidered on a torus. 

In the derivation of the continuum result, the starting point is the multi-particle phase 
space integration 


TV 

pae )=n / 

i=i 


dki 


2EM2nf 





( 20 ) 


where the energy of the Fth particle is E^, = y is an overlap factor can 

depend on the particle momenta. While the general result is Lorentz invariant, for simplicity 
we work in the rest frame, where P = 0. For energy just above threshold, E > Pxh, all 
of the individual particles are non-relativistic. This fact can be used to treat the overlap 
factors as constants, which is hne provided there is no orbital angular momentum. The iV-th 
particle’s momentum, say, is entirely hxed by three-momentum conservation. The energy 
conserving delta-function puts a constraint on the magnitude of the {N — l)-st particle’s 
momentum in the form 


dfcjv—1 5 


Pfc. — ~ J dk^-i Ftv-i d ^N-i — \/2Mf^_i [E — Pxh) 


( 21 ) 


To bound the N — 2 momentum integrals, we note that the maximum of any remaining 
particles’ energy is set by P — Pxh- Hence we have 


TV-2 


Pn{E) ~ -x/P — pTh [ 

i=l -^0 


d{ki)\/kf ^/E - E 


-3N-5 


Th 


( 22 ) 


which is the known scaling. 

For a torus, the starting point instead of Eq. 


IS 


TV .. oo 

iz»wi’=n73 E 


^ r 3 ;r(3) r 1^1 

_ r, / ~ Op* v/y 0 jp jp ^TV 

2Er 

i=l ni=—oo 


(23) 
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There can be additional suppression by integer powers of the available energy due to angular momentum 
considerations. For a state of orbital angular momentum £, we have an additional factor of {E — E^hY■ 
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where E is a non-integer index characterizing the energy of the A^-particle state. Notice 
in writing the energy-conserving Kronecker delta in Eq. ([23]), we have the total energy as 
the sum of the single particle energies. This is an approximation valid for weak interactions 
among the particles, because interactions modify the hnite volume spectrum of multi-particle 
states. The weight factor |Zjv(E)p is dimensionless because it is to be summed over the 
possible values for E, while Pn{E) has inverse mass dimension because it is to be integrated 
as a function of E. The case N = 2 appears to be somewhat special, because it can be 
evaluated more generally than at threshold. For N = 2, we have 


Z2m" 


1 

2ML3 


E 


1 

2Eh 


^E,M+E 



(24) 


where we treat one particle of mass m as light, and the other particle of mass M as heavy. 
We have not appealed to the threshold condition here. Because the overlap factor can only 
depend on n^, the sum collapses to a hxed value of determined by the size of E, and we 
hnd 

\Z^(E)\^<x - - . ( 25 ) 

AML'^ \/k‘^ + m? 

where the total energy is E ^ M -|- + m?. Unless there are momentum factors due 

to orbital angular momentum, no suppression is found near threshold. This occurs, for 
example, in the case of pion-nucleon coupled to an s-wave. For relative p-waves, our chiral 
perturbation theory result nicely lines-up with this general two-body analysis. In the low- 
energy expansion, oc - , due to the orbital angular momentum. 

Having handled the special two-body case, we now proceed to analyze the hnite volume 
phase space with generic particle number N >2. Considering specihcally near threshold, all 
particles are non-relativistic, and we choose the iV-th particle’s three momentum in Eq. fl2^ 
to vanish. Each sum over triplets of integers, n, can be changed to a sum over a radial 
integer n, with a multiplicity weight f{n) that counts the number of triplets n that can sum 
to a given n = y/^. For large n, we approach a continuum of states, and f{n) cx n. For 
small n, we take f{n) ~ 1 because corrections to this will give additional suppression at 
threshold. The remaining constraint due to the energy conserving Kronecker delta Exes the 
value of the (iV — l)-st particle’s energy. There is no change of variables factor here, however, 
we still have the multiplicity weight /(nAr_i), with un-i hxed by energy conservation. Thus 
we are led to 

f L \ N—2 VE—E^h ^ 2 

|Z^(F)|2~/ — V2M^_i(F-FTh) n E /(^.) oc - ^Th " , ( 26 ) 

\ ^ / i=l ni=0 

where we treat of generic size, but consider (F —FTh)T 1. From the intermediate 

expression, we can take the limit of large mode numbers and recover the correct scaling for 
continuum states. To achieve this, notice the behavior An/n —)■ dn, f{n) —)■ n, and we must 
accordingly take {E — FTh)T 3> 1 in the inhnite volume limit. For multi-particle states 
greater than two-body, we hnd the hnite volume result, Eq. fl26|) . maintains the vanishing 
of spectral weights near threshold, but not as sharply as for continuum states. 

For the case of a nucleon with pion interactions constrained by chiral dynamics, we expect 
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the spectral weights to behave as 




9a 


E^PL^ 


N-l 



2(N+i-2) 




(27) 


with E — h^Th ~ 2 ^' from threshold, the dimensionless ratio \k\/E^ is replaced by 

some complicated fnnction of this ratio, E{\k\/Ej^). The pre-factors can still be nsed to 
argne snppression of mnlti-particle states in the general case, provided that fL > 1, and 
rnTT-h S> 1. Contamination from mnltiple pions in the nncleon two-point fnnction shonld 
thns be largely irrelevant provided the smearing radins of the interpolating operator is small 
compared to the pion Compton wavelength. 


III. AXIAL CURRENT THREE-POINT FUNCTION 


Having derived the pion-nncleon and pion-delta contribntions to the nncleon two-point 
fnnction, we now tnrn onr attention to the determination of the nncleon three-point fnnction 
of the axial-vector cnrrent. The case of three-point fnnctions is complicated by the lack of a 
spectral representation. As a result, couplings to various multi-particle states have signs de¬ 
termined by the underlying dynamics. Cancellations can exist between these contributions, 
and we investigate whether chiral dynamics in hnite volume drives the axial three-point 
function above or below the single-nucleon contribution. In inhnite volume, multi-particle 
couplings drive the axial three-point function above the nucleon contribution 20j, which, 
however, would lead to an overestimation of the axial charge, qa- The same is true for 
discrete multi-particle states in hnite volume. 


A. Axial Current Correlator 

Our interest lies in the three-point function formed from the iso-vector axial-vector current 
inserted between two single-nucleon states. To determine the axial charge of the nucleon, 
we project the initial and hnal states onto vanishing three-momentum. Writing the current 
as we have the un-amputated three-point function 

G5„(r,i) = ^(0|OA,(f,T)J+(y,()Op0,0)|0>. (28) 

ky 

This correlation function depends on the source-sink separation, r, as well as the current 
insertion time, t, which we have assumed satisfy the hierarchy t > t > 0. Ideally one works 
in the limit in which r ^ f 3> 0, so that single-nucleon states can be cleanly isolated from 
both the source and sink. 

Inserting complete sets of hadronic states into the three-point function, we have the 
general expression 

(29) 

m^n 

Notice for ease below, we have factored out an overall exponential involving the nucleon 
mass and an overall operator overlap factor \Zo\^. Thus all of the energies are differences 
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relative to Mat, and all of the amplitudes Zn are ratios relative to jZaj- The ground-state 
contribution from the nucleon is isolated at long Euclidean time separations 

G5^(r,t) jZaj^e-^-^(N(0)jJ+jN(0)} + ■ • • , (30) 

where exponentially suppressed contributions from excited states have been dropped. In a 
non-relativistic notation, the axial-current matrix element in the nucleon has the form 

(N(d)IJ*^IN(0)) = 2ffA«*S„n, ( 31 ) 

where m is a rest-frame spinor, and is the covariant spin operator. The constant is 
the axial charge of the nucleon in the standard normalization. 

Considering now the general transition matrix elements of the axial current, we can write 
them in the form 

(njJ^^jm) = (32) 

with gnm as real-valued couplings, and cpnm as phases. The spinors appearing above charac¬ 
terize the multi-particle states. These are, in turn, constrained by the nucleon interpolating 
operator. The only overlap factors, that are non-vanishing correspond to states \n) 
which have vanishing total three momentum, and total angular momentum of Such a 
state can be described by a non-relativistic Pauli spinor u that merely characterizes whether 
the total angular momentum is aligned or anti-aligned with respect to some axis. Hence we 
dehne the transition matrix elements using the common spinor product keeping in 

mind that any differing normalization factors for multi-particle states are thereby absorbed 
in the gnm couplings. 

On account of Hermiticiy of the axial-vector current and time-reversal invariance, we 
have symmetric couplings, gnm = Qmm and antisymmetric phases ipnm = —^mn- The lat¬ 
ter condition implies that the imaginary part of the correlator vanishes. Taking this into 
account, we have the general form of the axial-three point function 

= 2u^ S^u\Zo\^e~^^'^'^gnmCOS(pnm\Zn\\Zm\e~^"^'"~^''e~^'^^ (33) 

n^m 

To analyze excited-state contributions to this correlator, we amputate the external legs by 
forming the standard ratio of three-point to two-point functions 


= (34) 

G(r) 

Isolating the single-nucleon contribution, this ratio has the general form 

^5m(g t) = 2u^S^u [gA + G'a('7-, t)] . ( 35 ) 

Notice the ground-state to ground-state coupling ^fQo is merely the nucleon axial charge 
5 'oo = 9A- While the single nucleon probability is not unity, it cancels out of the 

ground-state contribution in the ratio of correlators. The time-dependent function, Ga{t, t), 
thus characterizes the excited-state contamination; and, on account of the general form for 
the correlator given in Eq. fl5^ . it appears as 

GA{r,t) =^( 5 'nn- 5 ' 00 )|^npe"®"^+ gnm cos (pnm\Zn\\Zm\e~^"^''~^^e~^'^K ( 36 ) 

n>0 n^m 
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To arrive at this form, we made the approximation of neglecting any terms with four or 
more powers of overlap factors \Zn\- Such terms will be higher order in the chiral expansion, 
and arise from expanding out excited-state contributions to the two-point function in the 
denominator. The leading excited-state contribution from the denominator, however, gives 
rise to the (— 5 foo)-piece present in the diagonal sum. This subtracted term is analogous to 
wave-function renormalization. 

In order to determine the multi-particle couplings, we find it efficacious to rearrange 
Eq. fl36|) slightly. With redefined couplings, we have the compact expression 

GA{r,t)= Qnm\Zn\\Zm\e-’^-^^-^^e-^-\ (37) 

n,m—{0,0} 

for which the relation between couplings is 

Qnm ^nm{,9nn 9 a) T (1 ^nm)9nm ^nm- (2®) 


B. Chiral Computation 


One can determine the excited-state contamination in hnite volume using chiral dynam¬ 
ics. Previously the leading-order expression for the contamination function ^^(r, t) was 


obtained for the nucleon operator Ny in inhnite volume j2^. See Fig. 3 of that work for the 


corresponding Feynman diagrams. After some rewriting, the terms appearing in that result 
can be grouped according to the time dependence of the general expression given above in 
Eq. fl57)) . To match onto this expression, we must additionally convert from the continuum 
of states used previously to the discrete states permitted on a spatial torus. Finally we must 
compute pion tadpole diagrams that result from the effective theory representation of the 
(|, 0) 0 (0, i) lattice interpolating field in Eq. ([9]). These diagrams account for the fact that 
the effective theory nucleon operator N is perturbatively close to Ny. The required diagrams 
differ from those appearing in Fig. [1] merely by the insertion of an axial-current interaction 
along the nucleon line. It is a straightforward exercise to see that such axial-current opera¬ 
tor tadpoles cancel against the tadpoles in Fig. [T]when computing the ratio of three-point 
to two-point functions. This is not surprising, as the latter tadpoles contribute only to 
the wave-function renormalization. It thus remains to dissect the one-loop computation to 
obtain the multi-particle couplings. 


1. Nucleon to Excited State Matrix Elements 

The axial-current correlation function receives contributions from intermediate nucleon- 
to-excited state transitions. Because the overall nucleon mass dependence has been factored 
out of the three-point to two-point function ratio, these nucleon transition matrix elements 
are easily identified, because they lack dependence on either the source to current insertion 
separation, t, or current insertion to sink separation, r — t. First, we handle the matrix 
elements where the excited state is a pion-nucleon state. In inhnite volume, the contribution 
from nucleon to pion-nucleon matrix elements can be written simply in terms of the residual 
spectral functions. We have 

G0''(t,() = £dE [e-'"'"*' + ■ 

(39) 
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This expression can be reverted to an integral over the pion’s three momentum and this 
integral replaced by a momentum mode sum, as in Eq. ffTSj) . Carrying out these steps 
results in 


G 


N,nN 

A 



165 'aV|„ I i^aA 




+ e , 


(40) 


in finite volume. The weight factors \Zfi\, and \Z^\ have been determined from the nucleon 
two-point function calculation, and appear in Eq. fllOp . Comparing with Eq. fl37p . we see the 
first exponential above corresponds to contributions from N-to-nN transition matrix ele¬ 
ments, whereas the second are symmetrical vrA^-to-iV matrix elements. Given our definition 
of multi-particle couplings, we find that 


QN,Tr(n)N 



IS^AV I ryA 

27\/2 ” 


(41) 


where we necessarily label the pion-nucleon state by the pion momentum mode numbers. 

In infinite volume, intermediate state nucleon to pion-delta matrix elements can also be 
readily identified by their time dependence. We have the contributions 




(c t) 


(^5 'a + ^5'aa^ J dEp^^{E)[e + e 

£ dE^p^^(E)p,^(E + A) [e-<®+^)(-« + , 


which, on account of Eq. flTHl) , can be cast in the form 


(42) 


G 


N,itA 




25 

9a + ^ 5 'aa 


^A|2 


I ^ I 

27x/2' 


:3 — {En+^){E — t) _j_ g —(Es+A.)tj 

(43) 


in finite volume. Comparing with Eq. fl37p . we see the first exponential above corresponds to 
contributions from iV-to-vrA transition matrix elements, whereas the second contributions 
are symmetrical vrA-to-A matrix elements. Given our definition of multi-particle couplings, 
we find that 

_ f , \ IvAi I ,7 I f A A'\ 

QN,n{n)A — + ^S'AA I \Zf^ I 27 £2 

where pion-delta states have been labeled by the pion momentum mode numbers. This 
exhausts all nucleon to excited state matrix elements appearing in the axial-current corre¬ 
lation function. Before continuing, we remark that all multi-particle contributions to the 
three-point function above survive when either t S> 0 or r 3> t, but not both. 


2. Excited State to Excited State Matrix Elements 

The remaining contributions to the three-point correlation function arise from interme¬ 
diate excited state to excited state matrix elements of the axial current. There are contri¬ 
butions from both irN-to-nN matrix elements, as well as vrA-to-vrA matrix elements. In 
infinite volume, these such contributions appear in the expression 

£ dEp,^(E)e-’^- - + £ dEp,^(E) .[Af,) 
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FIG. 3. Excited-state contamination in the axial-vector three-point function of the nucleon, 
Eq. (j5ip . The contamination function t) is plotted as a function of the current insertion time 

t (measured relative to the midpoint between source and sink), for three values of the source-sink 
separation, r/a = 10, 12, and 14. Excited states with pions tend to drive the axial correlator 
upwards, which could potentially lead to an overestimation of qa of a few percent at most. 


Using Eq. ffT^ . the corresponding expression on a spatial torus appears in the form 




25 


gA\Z^\^e-^^^ + gA + -9 aa 


81- 


from which we can deduce the two couplings 


QttN^ttN and OttAjTtA I 9a T q.. 9aa 
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25 


81- 


( 46 ) 


( 47 ) 


These couplings are momentum independent. 

Finally, the one-loop infinite volume result contains TriV-to-vrA couplings. These appear 
in the expression 

t) = 1°° dE^J^AE)KAE + A) [e-'Jl-Oe-CE+AX + ^ 

( 48 ) 

The corresponding result on a spatial torus is immediate utilizing Eq. 

t) = J 2 \^n\\Zi\ ^e-^niT-t)^-iE^+A)t ^ g-(E-+A)(r-qg-E,,ij ^ ^49^ 

27\2 


and allows us to identify the final coupling, which is also momentum independent 


07rA7,7rA 
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27\/2 


gAN- 


(50) 


C. Axial Correlator Results 

We have determined the excited-state contributions to the three-point function of the 
axial-vector current in the nucleon, and thereby found the corresponding multi-particle cou¬ 
plings. Using the expressions presented above, we can write the excited-state contamination 
in the form 

A[T,t) = G^ +G^ +Gj^’ +G^ +Ga , ( 51 ) 
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where the ordering of terms is based on their expected size, and we have suppressed the 
dependence on source-sink separation, r, and current insertion time, t, on the right-hand 
side of the equation. In Fig. [3l we investigate the excited-contamination determined from 
chiral dynamics. Parameter values are taken to be the same as we used above for the two- 
point function: lattice parameters are a = 0.1 fm, and L = 48 a; the pion mass is the physical 
value, and the nucleon-delta mass splitting is A = 0.29 GeV, axial couplings are given values 
qa = 1.25, qan = 1-5, and = —2.25. While the hnal coupling has large uncertainties, 
our results are insensitive to the value due to exponential suppression. Results shown in the 
hgure are quite similar to those obtained previously in inhnite volume jl^. Contributions 
from intermediate states with pions tend to drive the axial correlator upwards, which would 
lead to an overestimation of the nucleon axial charge if such contributions are not accounted 
for. This hnding makes the consistent lattice underestimation of the axial charge more 
mysterious. 


IV. SUMMARY 


We consider the effect of pion-nucleon and pion-delta states on two- and three-point 
lattice QCD correlation functions. From the perspective of chiral dynamics, lattice QCD 
interpolating operators for the nucleon have dehnite transformation properties under chiral 
rotations, see [2^. As a result, these operators can be systematically mapped into effective 


field theory counterparts j25|. At leading order in the low-energy expansion, chiral correc¬ 


tions to the correlators appear with coefficients that are fixed due to the pattern of symmetry 
breaking j^. For lattice nucleon interpolating operators transforming as (|, 0) © (0, |), the 
interpolating field in the effective theory is not only perturbatively close to the heavy-nucleon 
held of chiral perturbation theory, the difference generates only tadpole corrections that 
do not affect multi-particle couplings. This justihes the considerations of an earlier study in 


inhnite volume [2^, provided the quark-level smearing of lattice interpolating operators for 


the nucleon is small on the scale of the pion Compton wavelength. 

For the nucleon two-point function in hnite volume, the derived couplings to intermediate- 
state pions above are shown to be in accordance with phase space available on a torus. While 
pion-nucleon and pion-delta spectral weights in hnite volume vanish less rapidly at thresh¬ 
old compared to inhnite volume, pion-nucleon contamination in the two-point function is 
estimated at the few percent level (a similar estimate is obtained in 2^), while pion-delta 
contamination cannot be discerned. For the nucleon three-point function of the axial-vector 
current, we similarly investigate the ehect of pion-nucleon and pion-delta contamination at 
zero momentum transfer. This three-point correlation function enables the lattice QCD 
determination of the nucleon axial charge, qa- As in inhnite volume, chiral corrections to 
the time-dependence of the axial correlator drive it upwards, with the dominant corrections 
arising from nucleon to pion-nucleon transitions. Thus for insufficient time between the 
source and current insertion, as well as between the source and sink, chiral corrections could 
potentially lead to a few-percent overestimation of the axial charge. This trend is oppo¬ 
site that encountered in lattice QCD calculations, and our study shows that the behavior 
is not likely due to pion-nucleon contamination. It would be interesting to study chiral 
contamination in other quantities determined from three-point functions, such as the quark 
momentum fraction in the nucleon, as results could help lattice practitioners better isolate 
the single-nucleon contribution. 
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